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ABSTRACT 

We study the effect of the magnetic field geometry on the oscillation spectra of strongly mag- 
netized stars. We construct a configuration of magnetic field where a toroidal component is 
added to the standard poloidal one. We consider a star with a type I superconductor core so 
that both components of the magnetic field are expelled from the core and confined in the 
crust. Our results show that the toroidal contribution does not influence significantly the tor- 
sional oscillations of the crust. On the contrary, the confinement of the magnetic field in the 
crust drastically affects on the torsional oscillation spectrum. Comparison with estimations 
for the magnetic field strength, from observations, exclude the possibility that magnetars will 
have a magnetic field solely confined in the crust i.e. our results suggest that the magnetic 
field in whatever geometry has to permeate the whole star. 
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1 INTRODUCTION 

It is known that magnetic fields play important role in the as- 
trophysical phenomena, such as supernovae, gamma-ray bursts, 
galaxy jet, and so on. Recently, there is a growing consensus in 
' explaining Soft Gamma Repe aters (SGRs) via the magnetar model 
, jPuncan & Thompsonlll993) . Magnetars are believed to be neu- 
' tron stars with strong magnetic field which is responsible for the 
" observed flare activity. Three giant flares, SGR 0526-66, SGR 
1900+14, and SGR 1806-20, have been detected so far. The peak 
luminosities are in the ranee of ~ 10'"' - 10'"' erg/s. This huge 
amounts of energy can be explained by the presence of a strong 
magnetic field whose stren gth is estimated to b e larger than 4x lO^'* 
Gauss for SGR 1900+14 faurlev et al.lll999h and of the order of 
8 X 10" Gauss for SGR 1806-20 jKouveliotou et al]|l998D . In ad- 
dition to the initial short and hard peak, these three events show a 
decaying tail lasting hundreds of seconds. In this late part careful 
analysis revealed the e xistence of characteristic Qu asi Periodic Os- 
cillations (QPOs); see IStrohmaver & Watts! ^Om for references. 
The oscillation frequencies in these QPOs are in the range of a 
few tenths of Hz to kHz and thought to be associated with crust 
torsional modes of magnetars. These signals might be the first ev- 
idence of a direct detection of crust oscillations of neutron stars: 
a first attempt to explain these frequencies by the crust torsional 
oscillations has been done b y using a model with dipo l e mag netic 
fields and realistic E oS (see lSamuelsson & Anderssonl ( 120070 and 
ISotani et all ( l2007al) l and this attempt was partially successful. On 



the other hand iGlampedakis et all ( l2006h suggested an alternative 
scenario to e xplain the ob servational frequencies via Alfven oscil- 
lations while lLeviiil ( l2006l) was the first to realize that the Alf ven os- 
cillatio ns of magnetars could be a continuum. More recentlv fLevinl 
( |2007|) pointed out that the edges or turning p oint of the continuum 
can be long-lived QPOs. In a recent paper, ISotani et al.l ( l2007bl) 
showed that the Alfven oscillations of magnetars could be a contin- 
uum which can explain the lower ones of the observed frequencies. 

However, it is quite difficult to construct realistic models for 
neutron star magnetic fields and/or to understand the details of their 
dynamics. In fact, the observational data suggest that the surface 
temperature of neutron stars with a strong magnetic field is not 
uniform (e.g.. lShowDe et alJ ( l2005l) : lHaberl et MilOO& i). This ob- 
servational evidence can be explained, for example, by the pres- 
ence of magnetic fields in the crust region with a strong toroidal 
components or a more complicated structu re than the poloidal one 
dGeppert et alj2006l :|P ons & Gepperll2007l) . Such crustal magnetic 
fields can be generated if the core is a type I superconductor. In 
this case, the magnetic fields are expelled from the core due to the 
Meissner effect after about 10* yearsQ. 

Recently, the analysis of the spectrum of timing noise for 
SGR 1806-20 and SGR 1900+14 has support ed the idea that th e 
core region can be a type II superconductor lArras et al]|2004l) . 
For the type II superconductor, even if an initial uniform magnetic 
field is present in the core, it will be split in thin magnetic fluxes 
tubes (fluxoids), which are parallel to the initial magnetic field. 
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These fluxoids interact with the neutron vortices and are driven by 
these vortices out side the core, tovy ards the crust on a timescale of 
10* — 10^ years / Rudermai] 1 2005h . In any case, it is too difficult 
to discriminate between all these models by using only the past 
observational data. This is due to uncertainties associated with the 
chemical composition of the interior and the crust of neutron stars 
as well as about the interior magnetic field configuration and its 
influence on the properties of neutron star. 

The aim of this paper is to calculate the axisymmetric crust 
torsional modes of magnetars with both poloidal and toroidal mag- 
netic field, where both the components are confined in the crust. 
Then our numerical results show that this magnetic configuration 
cannot explain the actual observational data of SGRs. In this paper, 
we adopt the unit of c = G = 1, where c and G denote the speed 
of light and the gravitational constant, respectively, and the metric 
signature is (— , +, +, +). 



2 EQUILIBRIUM CONFIGURATION 

The deformation from spherical symmetry due to the magnetic 
pressure would be small, because the magnetic energy is consid- 
erably smaller than the gravitational energy. So in this paper we 
neglect the effect of deformations induced by the magnetic field 
as well as the rotational deformations since the known magnetars 
rotate slowly. For these reasons, the equilibrium model of neutron 
stars can be described by a spherically symmetric solutions of the 
Tolman-Oppenheimer-Volkov (TOV) equations. The line element 
is given then by 



ds 
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and the 4- velocity is expressed as u'^ = (e"*, 0,0,0). In the fol- 
lowing, we assume an ideal MHD approximation and we consider 
an axisymmetric magnetic field generated by a 4-current J''. In or- 
der to determine the distribution of magnetic fields we need two 
basic relations: the Maxwell equations and the equation of motions 
that can be obtained by projecting the conservation of the energy- 
momentum tens or T'^'^ on to the hy persurface normal to u'^ (see 



momentum tens or 1'^ on to the hy persurtace normal to u'^ (s 
equation (11) in lSotani et al. I( l2007ij) ). These two equations are: 



pill' 



(2) 
(3) 



where is the Faraday tensor linked to the electromagnetic 4- 
potential by F^v = — A^,^. With the appropriate gauge 
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With and J* obtained from the equation l[2}, the equation ([S} 
for i-i = <j) can be written as 



Ar.ffA. 



A,rA^,9 = 0, (5) 

where Ar is defined as Ar = e'^'~^Ar,e sin 6. Thus Ar depends 
on A^ and we can set Ar — (^A^ (where ^ is a constant, see below 
for its meaning) since Ar is of the same order as A,f, with respect 
to the magnetic field strength. If we expand A^ as 

A^ = atj^j (r) sin 9dePt^, (cos 9), (6) 

then, the 4-potential A^ has the form 

A^ = (0, (e'^'^'^aej^jPtj^j , 0, af„ sin edgPt^,). (7) 
In the following, we will consider only dipole fields i.e. Im = 1- 



In order to obtain an expression for the 4-current , we can 
get additional equations from the equation ^ by setting fj, = r,d: 



X,r = Aifi^rJif, and x,e = A^^J^f,, 
where x is a function of r and 9 and 
1 
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Note that Cs is the sound speed, while in order to derive the above 
equations we have used the first law of thermodynamics, i.e, de — 
{e + p)dn/n, where n is the number density of particles. Using the 
integrability condition x,r9 ~ X,er together with equation ([S}, we 
obtain A^^rJ4,,e — A^^eJ^^r ~ 0. Thus we can see that is also 
function of A^ and it is possible to rewrite as = — fo — fiA^, 
where /o and /i are some constants. Comparing this last expression 
with equation ([9]l, we get the following expression for J^: 



J4, = — e 
47r 
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(11) 



For simplicity we neglect the contribution of /i 0. 

Setting /i = (/> in Maxwell equations ^ and by making use of 
equation dl lb we derive the following equation for ai 
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where ji = /or^ (e + p). For the exterior region, if we assume only 
dipole magnetic fields, i.e., C = 0, there is an analytic solution for 
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where fit is the magnetic dipole moment observed at infinity. The 
interior solution for ai should be determined by solving numeri- 
cally equation l ll2t and requiring that ai and ai' are continuous 
across the stellar surface. Hence, the components of the vector 
= B^/\/4iv, are given by 



Hr 



cos 9, 
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^From these expressions, we can see that the parameter is a con- 
stant associated only with the toroidal component of the magnetic 
field. Additionally the value of is chosen in the range where the 
function of ai has no node inside the star. Therefore there exists 
a maximum value of ^ named as (^max- Note that we can produce 
magnetar models with larger value of than (max- But for these 
models, the structure of the magnetic field becomes more complex 
since the sign of current is not constant in the star, and this could 
lead to possible instability of the magnetic fields. In other words, 
when the sign of current changes, even if the star remains stable 
the magnetic field structure might undergo a catastrophic instability 

^ Actually the term of /i can contribute in the equation I12t as the term of 
order £m = 3. However we use /i = since the contributions from terms 
with £f,i > 2 are negligibl e. Additionally, the omission of the /i term leads 
to the condition adopted in lKonno et alj j 19991) in the limit of <J ^ 0. 
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Figure 1. The three components, = R'^Hr/^it (for S = 0), Hg = 
B? He /fib (for e = 7r/2) and = B?H^/fjn, (for 9 = -k 12) of the 
magnetic field for a stellar model with 1.4Mq. The different lines corre- 
spond to different values of the constant f with C,max = 1.9. Note also that 
Hg is scaled by a factor ten. 



which will restructure the field. So we examine only the magnetar 
models with ^ ^ Cmax- 

In this paper we consider a stellar model with a superfluid core 
(type I superconductor) so that the magnetic field is expelled from 
the core by the Meissner effect. As a consequence, the magnetic 
fields are confined only in the crust region. In order to construct the 
model, we have to impose the condition that the magnetic fields are 
vanishing inside the core, i.e., we require that a\ should be zero at 
the basis of crust. This condition can be written as a\ ~ a.c.(r — 
Re), where Qc is a constant and Rc is the radius at the basis of the 
crust. Note that at r — Rc both Hr and vanish while Hg is not 
zero. In Figure[T] we plot the components of magnetic field, where 
Hr, He, and H^ are defined as Hr = R?Hr/ fJ-b, Hg = R^Hg //ib, 
and H^ = R?H^/ ^b- The actual details of the used stellar models 
will be described in the next section. It is worth noticing, from these 
figures, that the 9 components of magnetic fields are not smooth at 
the stellar surface, because we assume an exterior magnetic field 
with only a poloidal component, which means C = for r ^ R 
i.e. no toroidal component outside the star. This choice produces 
a current on the star surface. The presence of a surface current is 
found also in many simulations dealing with the evolution and th e 
stability of the magnetic field (see Braithwaite & Spruit, l l2006h l. 
Howevere we do not consider the existence of the surface current 
in this paper, since the surface current dose not affect the behavior 
of the torsional oscillations within our analysis. 



3 NUMERICAL RESULTS 

In this paper, we adopted polytropic equations of state (EoS) de- 
scribed by the following relations 

(71 \ ^ P 
— and p = nmb + — -, (17) 
no/ r — 1 

where — 1.66 x 10"^'* g and no — 0.1 fm~^, and we have 
set r = 2.46 and K = 000936 by fitting the tabulated data for 
EoS A IPandharipandd l 197ll) . With this equation of state and with 
the fixed value of the density at r = R^, p = 2.4 x lO^** g/cm'^, 
the maximum mass is Mmax — 1.654M0 . Here we study two rep- 
resentative stellar models; one with M = 1.4M0 and the other 
is the one with the maximum allowed mass for this EoS i.e. the 
one with M = 1.654i\f©. For these two stellar models, the radius, 
the relative crust thickness, and the compactness are R = 10.35 
km, Ar/R = 7.34 %, and M/R = 0.200 for the model with 
M = 1.4Mo while the model with M = 1.654M0 has i? = 8.95 
km, Ar/R = 4.15 %, and M/R = 0.273. Finally, for simplicity 



the shear modulus in the crust, /i, is assumed to be given by a rela- 
tion of the formula fi = v^p, where Vs is the speed of shear wave s 
with a typical value Vs ~ 10* cm/s ('Schumaker & Thomel |l983h . 
In general the frequencies of torsional oscillations depend on the 
EoS in the crust r egion since the shear modulus also depend on the 
EoS. However, in lSotani et alj l l2007ah we have used two different 
realistic EoS for the crust and our results show that the freuqen- 
cies of fundamental torsional oscillations are almost independent 
of them. Anyway in this paper, as first approximation, we use the 
above simple formula for shear modulus. 

We study axial perturbations in the Cowling approximation of 
the previous stellar models. We focus only on these type of oscilla- 
tions of the crust (torsional) since they do not induce density per- 
turbations, which means that significantly lower energy is needed 
for their excitation. On the other hand the absence of density per- 
turbations leads to minimal metric perturbations which justifies our 
choice in using the Cowling approximation. Under these assump- 
tions the pe r turbatio n equations are reduced to equation (51) in 
ISotani et al. I ( l2007ij) . By studying this equation one can observe 
that, if the background magnetic fields are axisymmetric, the ax- 
isymmetric torsional oscillations are independent from the compo- 
nent H^ of the magnetic field. Thus we can calculate the eigen- 
frequencies by usin g the perturbation equations (69) and (70) in 
ISotani et al. I j2007ah with exactly the same boundary conditions, 
such as y^r = at the basis of crust (r = Rc) and at the stellar 
surface (r = R), where y{r) describes the angular displacement 
of the stellar material. How ever, we have to replace ji in equation 
(70) in lSotani et al.l ( l2007i^ by h = ji +C^e"2*ai/(47r), because 
for the elimination of the term of ai" we have used equation J12b . 

The numerical results are shown in Figure (2] where the fun- 
damental 1 = 2 torsional modes are plotted as function of the nor- 
malized magnetic field B/B^^ (B^ = 4 x 10^^ Gauss). ^From 
these figures, it is clear that the frequencies depend weakly on the 
presence of the toroidal magnetic fields, even if the strength of the 
toroidal field increases i.e. the value of ( becomes larger. This weak 
dependence can be understood if we recall that the magnetic field 
lines of the toroidal component are parallel to the direction of fluid 
motion. But the most important result is that, although it seems 
that qualitatively the effect of the magnetic field is similar to the 
results in Sotani et al. (2007a), there is a significant quantitative 
difference: the influence of the magnetic field on the axisymmet- 
ric torsional frequencies becomes apparent for considerably lower 
strengths of the magn etic field, i.e. already for B /B^ ~ 0.1. No- 
tice that the results in lSotani et alj ( 12007 al) show that the effect of 
the magnetic field can be seen for B/B^ > 1. It is natural that 
when the magnetic field is confined in the crust, magnetic field be- 
comes considerably stronger than the one that permeates the whole 
star for the same strength of the exterior magnetic field. Actually, 
since the magnetic field is confined in the crust which has thickness 
Ar ~ 0.03 - 0.12i? ^ 0.5 - 1.5km, the density of magnetic field 
lines becomes higher. Also if we compare the results for the two 
stellar models under discussion we observe that thinner is the crust, 
then earlier the effect of the magnetic field becomes pronounced. 
For example, for the model with M = 1.654M0 (Ar/R = 4.15 
%) the effect of magnetic field is apparent for B/B^i ~ 0.07. For 
the other model (M = 1.4M0 and Ar/R = 7.34 %) the effect of 
the magnetic field becomes pronounced later, i.e., B/B^ ~ 0.1. 

The above results are derived by using the same procedure as 
in lSotani et all j2007ah , i.e. we decompose the system into spheri- 
cal harmonics and truncate the (£ ± 2) terms. This truncation could 
introduce small errors in the results, thus in order to check the va- 
lidity of our results we have used an alternative technique based 
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Figure 2. With M = IAMq and 1.654M0, the frequencies of the 
fundamental torsional modes of £ = 2 as functions of B/B^i, where 
B/j, = 4 X 10^^ Gauss. The different lines correspond to different val- 
ues of (, where (max ~ 1-9 for M = 1.4Mq and (max ~ 3.0 for 
M = 1.654M0. 
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Figure 3. The torsional mode frequencies as functions of the magnetic field 
strength for ( = 0. The continuous line shows the results from the 1-D trun- 
cated mode calculation on which we omitted the coupling to £±2 equations 
while the thick dots correspond to results from 2D numerical evolutions of 
the perturbation equations. In this figure two vertical lines has been drawn 
corresponding to the estimated magnetic field strengths for the two SGRs 
that is B Ki 4 X 10^'' Gauss for SGR 1900-H14 and B Ri 8 X 10^'' Gauss 
forSGR 1806-20. 



on 2D numerical evolutions of the perturbation equation describing 
the dynamics of the oscillatin g crust. The 2D m ethod used here is 
similar to the one described in lSotani et d]( l20b7b.) and the evolu- 
tion equations were structurally similar although the various back- 
ground terms have been modified to include terms related to the 
shear modulus, fi. The 2D evolution was constrained only in the 
crust and in the first quadrant i.e. ^ 9 ^ ir /2 and the following 
boundary conditions have been imposed: (a) y^r = at r = i?c 
and R, (b) y,0 ^ at 9 = 0, and (c) 3^ = at 61 = 7r/2. The 
results of the two approaches for — are plotted in in Figure 
[3] The continuous line corresponds to results derived by using the 
truncated equations and treating the problem as a boundary value 
one while the thick dots correspond to results produced by the 2D 
numerical evolution. The agreement between the two approaches is 
apparent, which suggests that for this type of problems the trunca- 
tion of the (£ ± 2) equations does not affect at all the results. 

Furthermore, we performed similar calculations for stellar 
models with different realistic EoS, but the behavior of frequen- 
cies, as function of a magnetic field, is similar to that for the EoS A 
that we use here. 



4 CONSTRAINTS ON THE MAGNETIC FIELD 
GEOMETRY 

We have studied the torsional oscillations for the magnetic field 
distribution where the magnetic fields has been "expelled" from the 
superfluid core and confined in the crust. We have also assumed that 
the magnetic field has both poloidal and toroidal components. The 
study shows that the presence of the toroidal field affects marginally 
the oscillation frequencies of the torsional modes, which can be 
explained because the perturbations propagate along the toroidal 
field lines. The toroidal field might play a more significant role in 
the case of non-axisymmetric perturbations or if we include the 
coupling to polar perturbations. 

The main results of this work is summarized in Figure[3] Our 
calculations show that the effect of magnetic field on the torsional 
modes, when it is confined only in the crust, can be observed for 
magnetic field strengths below 10^^ Gauss. Moreover, this obser- 
vation has a significant effect on understanding the geometry of the 
magnetic field of magnetars. The reason is the following. With this 
magnetic field configuration there are no additional oscillational 
frequencies such as Alfven oscillations (see Sotani et aL t2007h)), 
because the magnetic field does not exist in the core. Thus the 
1 = 2 torsional oscillation is lowest oscillational frequency. On the 
other hand the estimated magnetic field strengt hs are of the order 
of a f ew times 10^'' Gauss (Hu rlev et alJl 19991 : iKouveliotou et al.l 
Il998l) or even higher ( iNakagawa et al. 2007). This suggest that if 
the magnetic field was constrained to the crust the lowest possible 
frequency for the signals from SGR 1806-20 is about 40 Hz. Then, 
this magnetic field distribution can not explain the observed fre- 
quencies 18, 26 and 29 Hz. So this type of magnetic field geometry 
cannot be present in the observed magnetars and most probably it 
has to be excluded in general for the case of strongly magnetized 
neutron stars. Alternative scenarios which might include coupling 
with the polar modes should be excluded since crustal perturba- 
tions of polar type are typically of conside rably higher frequencies 
( IStrohmaverlll99ll : iWoulidis et alj|2007l) and cannot explain all 
three frequencies. 
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APPENDIX A: POLOIDAL AND TOROIDAL FIELDS 
PERMEATING THE STAR 

In this appendix we study the effect on the torsional mode spectrum 
of the toroidal comp onent added on the p oloidal magnetic field con- 
figuration studied in lSotani et alj ( 12007 ah . In this case the magnetic 
field permeates the whole star but we study only the oscillations of 
the crustal part of the star. Although this is somehow a truncation 
of the full problem, it may be a good approximation in studying 
the effect of the toroidal component (se e the next paragraph for 
an argument in favour of this). Following lSotani et al.l ( l2007al) . the 
equilibrium configuration of the magnetic field is constructed by 
imposing the regularity condition, ai ~ Qc^^, at the stellar center, 
and by requiring the continuity of ai with the analytic exterior so- 
lution l |13t . The contribution of the toroidal field is set to zero for 
r ^ R. The components of the magnetic field are plotted in Figure 
lAll in a similar fashion as in Figure[T] 

The torsional oscilla tion frequencies are calculated following 
the same procedure as in ISotani et al.1 ( l2007ij) , i.e., 1-D truncated 
mode calculation using the same boundary conditions, i.e., y,r = 
air — Rc and r — R. Actually, the boundary condition atr = Rc 
is not quite correct because in this wa y we ignore the coupling 
with the fluid core jSotani et al.ll2007d) . Thus the proper bound- 



ary condition is Hry^r ' = [/i + H^Hr] ^.'v^'', which can be 



(+) 



B/B^. iFvom 



written as Vy'/y;^' ~ P /{P + 1), where (3 ■ 
this relation it is evident that for magnetic fields with strength 
around B/B^ ~ 0.1 — 0.3 the coupling with the fluid core is 
not so important for the study of pure crust torsional oscillations. 
The corrections to our approximate boundary condition arise only 
for B/Bfi > 1. Thus approximately we adopt the simple bound- 
ary condition y^r = at the basis of crust, which will be used to 
demonstrate the effect of the toroidal component on the spectrum 
of torsional oscillations. 

In Figure IA2I we show the frequencies of the fundamental 
1 = 2 torsional modes as function of the normalized mag netic 
field strengh, B/B^. As mentioned in ISotani et alj j2007ah . the 
effect of magnetic field on torsional modes becomes important 
for B ~ B^. This is expected since for B > the Alfven 




Figure Al. The three components, Hr = R'^Hr/^ib (for 9 = 0), 
Hg = R'^ He /fit (for 9 = 7r/2) and H^, = R^H^/^n, (for 6» = 7r/2) 
of the magnetic field for a stellar model with IAMq. The different lines 
correspond to different values of the constant C, with C,max = 0.31. 
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Figure A2. The frequencies of the fundamental torsional modes with 1 = 2 
as functions of B/B^ with M = \AMq and 1.654Mq. The different 
lines correspond to different values of where Crnaa; ^ 0.31 for M = 
1.4M0 and Cmaa: ~ 0.27 for M = 1.654M0. 



speed, ua = B /{A-npY^"^ , becomes larger than the shear speed, 
Ug = (/i/p)^''^. ^From these figures, it is clear that the presence 
of the toroidal magnetic field affects weakly the frequencies of the 
torsional oscillations. 

Concluding, the presence of a toroidal magnetic field compo- 
nent does not affect significantly the spectrum of the crust oscilla- 
tions neither in the case when the poloidal field permeates the star 
nor when it is confined in the crust. This suggests that although we 
are able to exclude the presence of crust confined magnetic fields 
in magnetars, most probably we will not be able via the QPOs ob- 
servations to identify the presence of a toroidal component in the 
magnetic field. 



